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ABSTRACT 

The galaxy Luminosity Function (LF) has been estimated in the rest frame B luminosity at < z < 
1.25 and at 1700 A for 2.5 < z < 4.5 from deep muhicolor surveys in the HDF-N, HDF-S, NTT-DF. 
The results have been compared with a recent version of galaxy formation models in the framework 
of hierarchical clustering in a flat Cold Dark Matter Universe with cosmological constant. The results 
show a general greement for z <; 1, although the model LF has a steeper average slope at the faint end; 
at z ~ 3 such feature results in an overprediction of the number of faint {Iab ~ 27) galaxies, while 
the agreement at the bright end becomes critically sensitive to the details of dust absorption at such 
redshifts. The discrepancies at the faint end show that a refined treatement of the physical processes 
involving smaller galaxies is to be pursued in the models, in terms of aggregation processes and/or stellar 
feedback heavily affecting the luminosity of the low luminosity objects. The implications of our results 
on the evolution of the cosmological star formation rate are discussed. 

Subject headings: galaxies: fundamental parameters - galaxies: evolution ~ galaxies: formation 



1. INTRODUCTION 

The evolution of the luminosity function (LF) is a fun- 
damental probe for cosmological theories of galaxy forma- 
tion. Indeed, the red and blue/UV luminosities are related 
to the galaxy stellar mass and star formation rate, respec- 
tively. Thus the evolution of the LF reflects the differential 
contribution of different galaxy types to the cosmic history 
of mass growth and of star formation which are the main 
outcomes of hierarchical models for galaxy formation. 

Adopting a Schechter fit to the LF, recent estimates 
(Zucca et al. 1997, Marzke & da Costa 1997, Folkes et 
al. 1999) point toward a relatively steep power law at the 
faint end (with a power index ~ 1-2), although an excess 
of dwarf blue galaxies relative to the Schechter fit is found 
at the faint end. 

In the intermediate redshift range (up to z ^ 1), first 
steps toward the evaluation of the LF evolution were un- 
dertaken by the Canada France Redshift Survey (CFRS, 
Lilly et al. 1995) and by the Autofib Redshift Survey (El- 
lis et al. 1996, Heyl et al. 1997). These surveys have 
shown some increase in the number density of the fainter 
population together with some increase in the luminosity 
of the brighter blue fraction. These two effects are respon- 
sible for the strong increase with z of the average cosmic 
UV luminosity density (by a factor 5-10) in the redshift 
interval < z < 1. 

Deep multicolor surveys of galaxies represent an effec- 
tive way to explore the galaxy distribution in the redshift 
interval 1 < z < 5. Successful spectroscopic confirmation 
was obtained for the brightest fraction by Stcidcl ct al. 
(1996,1999). For the bulk of the population reliable pho- 
tometric redshifts are currently used (e.g. Connolly et al. 
1997; Giallongo et al. 1998, Fernandez-Soto, Lanzetta & 
Yahil 1999; Fontana et al. 2000). 



The average UV luminosity density resulting from these 
studies shows a fiat distribution at z > 1 that extends up 
to z ~ 4.5. Such a result represents a strong constraint to 
hierarchical galaxy formation theories that predict a signif- 
icant contribution by a large population of small objects. 
In a previous paper we showed that when the effects of 
the sample magnitude limits are included in the models, 
a steadly decline in the predicted UV luminosity density 
appears in contrast with the observed behavior (Fontana 
et al. 1999). 

Attempts to clarify the contribution of different (in mass 
and SFR) galaxy populations to the global UV luminos- 
ity density, requires an evaluation of the high redshift LF 
down to the faintest accessible magnitudes. In this Let- 
ter we present a first estimate of the intermediate/high- 
z luminosity function using a composite deep multicolor 
sample of about 1200 galaxies with reliable photometric 
redshifts down to Iab = 27.5. The depth of this data sam- 
ple, together with its relatively large area, allows a direct 
comparison of the observed LF shape and evolution with 
theoretical predictions. Such comparison is performed us- 
ing a semianalytic implementation of recent hierarchical 
models of galaxy formation. 

2. THE DATA SAMPLE 

The analyzed dataset covers a wavelength range from 
the UV to the K bands, and observations were taken from 
ground based telescopes and from the HST as well. The 
first field, known as New Technology Telescope deep field 
(NTTDF) consists of an area of 4.84 arcmin^ where op- 
tical and near-IR UBVRIJK observations have been taken 
at the ESO NTT telescope with various instrumentation 
(SUSI, SUSL2, SOFI). Further details about the muhi- 
color catalog can be found in Arnouts et al. (1999) and in 
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Fontana et al. (2000). 

We have also used the Hubble Deep Fields North and 
South catalogs provided by Fernandez-Soto et al. (1999), 
with an overall area of 3.92 arcmin^ and 4.22 arcmin^ re- 
spectively. After appropriate selections to remove con- 
tamination by stars and low signal-to-noise ratio (S/N) re- 
gions, we applied our photometric redshift code to the data 
down to Iab = 27.5. A detailed description of this proce- 
dure, along with photometric z catalogs, can be found in 
Fontana et al. (2000). 

In addition, in the present evaluation of the luminosity 
function, we have also considered one of the main system- 
atic errors affecting the estimates of the galaxy total abso- 
lute magnitudes at various z, which is the surface bright- 
ness cosmological dimming effect. To correct for this effect 
one has to recover the total galaxy emitted flux. To esti- 
mate the systematic losses associated with any recovering 
procedure, we performed simulations of synthetic galaxies 
(having exponential intensity profiles with an ellipticity of 
0.5) with different apparent magnitudes, as seen with the 
appropriate S/N in the HDF images. Then we compared 
the total input flux with the one obtained by extrapolating 
the intensity profiles as computed in the SEtractor package 
(Bertin & Arnouts 1996). The simulations show that the 
difference A/ between the input and the measured appar- 
ent magnitudes increases from 0.1 to 0.25 when the input 
magnitude increases from 7 ~ 25 to / ~ 27.2. For these 
reasons, although the catalogs in the HDFs go deeper, we 
will confine the analysis of the LF at magnitudes / < 27.2 
in the HDFs (and 25.7 in the NTTDF) where errors in 
the estimate of the total magnitudes are small and do not 
affect our main results on the shape of the LF. 

3. ESTIMATING THE LUMINOSITY FUNCTION 

Several methods are available in the literature (see Ef- 
stathiou, Ellis & Peterson 1988 for a discussion): here we 
choose to adopt the classical l/Vmaa- estimator (Schmidt 
1968) jointly with the Sandage,Tanunami & Yahil (1979) 
maximum likelihood fit for a Schechter function. In the 
1/Knaa; method, for any given redshift bin (21,^2) an ef- 
fective maximum volinnc is assigned to each object. This 
volume is enclosed between zi and Zup, the latter being de- 
fined as the minimum between Z2 and the maximum red- 
shift at which the object could have been observed given 
the magnitude limit of the sample. 

Combining data from separate fields with different mag- 
nitude limits, we then compute the galaxy number density 
4>{Mb,z) in every (A^;, AMb) bin as follows: 
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where w(j) is the area in units of steradians correspond- 
ing to the field j and N is the total number of objects in 
the chosen bin. The number of fields involved in the sum 
over index j is restricted to the ones with a faint enough 
magnitude limit for the i-th object to be detected. In 
this way each galaxy has a different Zup{j) in every field 
and the overall volume available for this object is obtained 
summing the corresponding Vmax{j)- 

On the other hand the STY technique, once assumed 



a Schechter behaviour for the LF, is an endeavour to 
maximize the likelihood of representing the observed set 
of galaxies with the best fit parameters of the Schechter 
function. Assuming that, in an appropriately thin bin of 
redshift, the number of sources with redshift between z 
and z + dz, and with absolute magnitude between Mb 
and Mb + (IMb can be factorized as (t){M b , z)dzdM b = 
p{z)ip{MB)dzd Mb, where p{z) is the density of galaxies 
at redshift z, that can be considered constant within the 
bin. The function iP{Mb) is taken to have a Schechter 
form, parametrized by the characteristic absolute magni- 
tude M^ , by the logarithmic slope at the faint end a and 
by the normalization V'* • 

If we choose an appropriate redshift bin (21,^2), it is 
possible to give an estimation of a and Mg considering 
the probability density pij to find a galaxy with redshift 
in the range Zi,Zi + dz and absolute magnitude between 
Mi and Mj -|- dM in the j-th magnitude-limited field and 
maximizing the likelihood of observing the set of galaxies 
that comes from the surveys, wich is simply given by the 
product of the single probabilities: 



p{zi)ip{Mi)dzdM 

rM;;i(z), 



Nj Nj 

a;(j) Q piz)^dzjT--'" m)dM 

(2) 

uj{j) beeing the field area in steradians, and Mll-^{z) the 
absolute magnitude value that the i-th object, if detected 
with the magnitude limit m^^^ in the j-th survey, should 
have at that redshift. This value clearly depends on the 
details of the spectrum of each object. Here j runs over 
the number of fields where the i-th galaxy can be detected, 
each one containing Nj objects. 

The value of (f)* is then obtained summing the density 
in the z,M space for every galaxy, taking account of all 
the fields where it could be detected (i.e. the fields with 
enough bright magnitude limit): 



^(i) / -r^^ / V(M)dM 



jdetect 



As for an appropriate selection in magnitude of the sam- 
ple, it is important to bear in mind that once a rest frame 
wavelength A is chosen (e.g. the 4400 A band), this cor- 
responds to different observed wavelengths A(l + z) when 
z runs inside the redshift bin (zi, ^2)- Since our aim is to 
choose appropriately a complete subsample selected in the 
rest-frame, we ought to take into account the different red- 
shifts of each galaxy and the details of the spectra as well, 
to be sure that we are selecting objects in a coherent way. 
One major concern is fixing the bins in redshift in a suit- 
able way for matching the observed bands once the wave- 
length displacement is taken into account. Subsequentely 
the right selection criterion can be achieved by comparing 
the photometric limit with the A(l + z) magnitude value 
taken from the spectrum. This implies a limiting magni- 
tude < 27.2 at a wavelength of 4400(1 + z) kiov z< 1.25 
and < 27.2 at 1700(1 -|- 0) A for ^; > 2.5 in the HDFs. 

4. THE COSMOLOGICAL EVOLUTION OF THE LF 



Poli et al. 



3 



-1 



-4 



- n 



1 \ I \ r 



M,A 



1,0 Ho = 50 



0.2<z<0.5 




a = -l.l 
M,= -21.36 
0.=O.OO59 
N=188 




-4 



nM,A=l'0 H„ = 50 



- 0.5<z<0.75 




a = -l.l 

M. = -21.0 
0, = 0.0091 
N = 251 



-3 



-4 



«M,A-1.0 H, = 50 
0.75<z<1.25 




a = -1.24 

M, = -21.45 
0, = 0.0045 
N = 409 



-24 



-22 



-20 



18 



16 



24 



-22 



-20 



18 



16 



Fig. l.-Rest frame B luminosity function of field galaxies in different redshift bins and for two cosmological models. Left panel: The continuous 
curves are the Schechter LFs resulting from the Maximum likelihood fit to our composite galaxy sample. Dotted curves are the Schecter 
functions derived from the CFRS survey (Lilly et al. 1995), while the dashed curves are derived from the Autofib redshift survey (Heyl et al. 

1997). Right Panels: Continuous curves arc the CDM model predictions discussed in the text, including dust absorption with a SMC extinction 
curve (different extinction curves used in Fig. 2 do not produce appreciable changes in the B band LFs). Dashed curves are the CDM dust- free 
model predictions. 



4 



The luminosity function at the rest frame B magnitude 
in the AB system, Mg, is shown in Fig. 1 for three redshift 
bins (0.2-0.5,0.5-0.75,0.75-1.25). The left panels refer 
to a critical universe with Ojvf = 1, = and Hq = 50 
km s""'^ Mpc"-'^, while the right panels refer to a flat uni- 
verse dominated by the cosmological constant CIm = 0.3, 
Qa = 0.7 and Hq = 70 km Mpc"^ The rest frame B 
magnitude is computed from the best fit theoretical spec- 
tral energy distribution used to derive the photometric red- 
shift. The magnitude limit of the sample has been eval- 
uated at the same rest-frame wavelength centered in the 
rest frame B band, which roughly corresponds to the ob- 
served V,R,I bands for the three redshift bins, respectively. 
In addition we show in Fig. 2 the 1700A rest-frame LF in 
the bins 2.5 < z < 3.5 and 3.5 < ^; < 4.5 for the same 
cosmologies. 

The best fit values of the Schechter parameters obtained 
from the maximum likelihood method are summarized in 
table 1 together with the relevant sample parameters for 
both the considered cosmologies. In the highest redshift 
bin (3.5 < z < 4.5) the small number of objects prevented 
an accurate estimate of the Schechter parameters. The 
Schechter curves corresponding to the tabulated fitting pa- 
rameters are shown in Figs. 1,2 only for the case of the 
critical universe (left panels) for comparison with previous 
LFs computed in the same cosmology. In particular, the 
Schechter LFs from the Canada France Redshift Survey 
and the Autofib survey are also shown for z < 1.25 as dot- 
ted and dashed curves, respectively. An overall agreement 
with the spectroscopic data is evident at bright magni- 
tudes, supporting the reliability of our LF photometric 
estimation. The variance between spectroscopic surveys 
could be due to different k-corrections applied for the es- 
timate of the rest-frame blue magnitudes from samples 
selected in different bands. As discussed in Sawicki et al. 
(1997), the uncertainties in the photometric redshift esti- 
mates do not affect appreciably the behavior of the LF: 
suitable Monte Carlo simulations have shown that pertur- 
bations arising from redshift uncertanties result in effects 
smaller than one sigma error bars in the 1 /Vmax estimator 
and in small changes in the best-fit Schechter parameters. 
In our case the good match with CFRS leaves little space 
for this kind of substantial changes in the distribution. 

It can be seen that from z ~ 0.2 up to 2; ~ 1.25, there 
is no evidence of a significant trend with redshift in the 
faint-end slope parameter a, which remains close to its 
local value a ~ —1.2 (Zucca ct al. 1997). Also, the char- 
acteristic magnitude shows a mild brightening with the 
look-back time. At high redshifts a steepening effect is ev- 
ident in the 1700A LF (Fig. 2, left panel), where the slope 
parameter reaches a value of —1.37, consistent within un- 
certainties with the a = —1.6 found by Steidel et al. (1999) 
in an analogous redshift bin. We note however that in this 
high redshift bin the slope is weakly constrained by the 
present depth of the galaxy sample. 



5. A COMPARISON WITH HIERARCHICAL MODELS FOR 
GALAXY EVOLUTION 

We compare our data with our rendition of the semian- 
alytic models. The structure of such rendition is described 
in Poll et al. (1999). We updated our model to implement 
some of the most recent improvements inserted in the re- 



cent versions of scmianalytic models, following the lines of 
Cole et al. (2000). In particular we now adopt the Lacey 
& Cole (1993) dynamical friction timescale, the new star 
formation and feedback recipes and the new modelization 
for hot gas distribution implemented by the above authors. 
This allows us to obtain a TuUy-Fisher relation in reason- 
able agreement with observations together with a good fit 
to the local luminosity function, the two observables that 
have been used to calibrate the free parameters governing 
the star formation and the feedback processes. In addition 
we have included dust absorption, which is modelled as in 
Poll et al. (1999). We have checked that the predictions 
of our model agree with those from Cole et al. (2000) 
for both the local observables and the global cosmic star 
formation history. 



Table 1. 

Parameters of the Schechter function 



fits 



z range 


a" 




<j>* 


N 


02-0.5 


-1.18 ±0.05 


-21.36 ±0.37 


0.0059 


188 




-1.19 ±0.05 


-21.03 ±0.38 


0.0086 




0.5-0.75 


-1.18 ±0.06 


-21.00 ±0.29 


0.0091 


251 




-1.19 ±0.06 


-20.75 ±0.30 


0.01 




0.75-1.25 


-1.24 ±0.06 


-21.45 ±0.24 


0.0045 


409 




-1.25 ±0.05 


-21.38 ±0.25 


0.0041 




2.5-3.5 


-1.37 ±0.20 


-20.72 ±0.37 


0.0025 


187 




-1.37 ±0.19 


-20.84 ±0.37 


0.0023 





" The second row for each z bin refers to the Qm = 0.3, 
rt\ = 0.7 cosmology. 

M* is Mb (AB) except in the z = 2.5 - 3.5 bin where 
it is computed at 1700 A, M1700 

To compare with the data we focus on the flat ACDM 
(Cold Dark Matter) cosmology. Such choice is favoured by 
a significant amount of recent observational results, from 
the high-fraction of the baryons-to-DM ratio in galaxy 
clusters (see, e.g.. White & Fabian 1995) to the high- 
redshift SNe (see Perlmutter et al. 1999). When inserted 
into galaxy formation models, the above set of parameters 
also yields a cosmic star formation history in reasonable 
agreement with observations (see Fontana et al. 1999), in 
contrast, e.g., with the Standard CDM with Qm = 1 and 
Qa = 0. 

Fig.l and Fig. 2 (right panels) show a comparison of 
the observed luminosity functions derived from our com- 
posite sample and our CDM model. First we note a general 
agreement for Mb < —20 and z <^ 1. Both data and pre- 
dictions show weak luminosity evolution resulting from a 
balance between a mild decrease with z of the number of 
massive objects and the increase of their star formation 
activity and hence of their blue luminosity. At the same 
time, the normalization of the LF rises by a factor of ^ 2. 

At low luminosities the theoretical LF appears steeper 
than the observed one, with an excess at the faint end 
that becomes larger with increasing redshift. In addition, 
at z > 3.5 the model tends to slightly underestimate the 
magnitudes of bright galaxies when the model with larger 
dust extinction (~ 1 mag, in agreement with Pettini et 
al. 1998) is considered. Such results show that a refine- 
ment of the models is required, both in terms of dynamical 
processes (like merging of satellites in common halos, see 
Somerville, Primack, & Faber 2001) and in terms of stellar 
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Fig. 2.-Rest frame ultraviolet luminosity function at z ~ 3 and 2 ~ 4. Filled symbols are derived from the spectroscopic survey by Stei- 
del et al. (1999). Left Panel: Continuous curve is derived from the Schechter Maximum Likelihood fit as in Fig.l. Right Panel: Symbols as 
in the corresponding panel in Fig.l except for the dotted or long-dashed curves which refer to the Milky Way and Calzetti extinction curves, 
respectively. 
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processes (like the feedback from Supernovae, appreciably 
affecting the LF at the faint-end). Indeed, the apparent 
agreement of the predicted vs. observed UV luminosity 
density at z > 3 (Fontana et al. 1999) results from the bal- 
ance in the LF between the excess of the predicted dwarfs 
and the deficit of predicted bright galaxies, as shown in 
the right panels of fig. 2. 

This shows that the 2:-resolved LFs constitute a particu- 
lary powerful way of constraining the two most uncertain 
processes in the theoretical modelling, i.e., the dust ab- 
sorption affecting the prediction of the number of bright 
galaxies and the Supernovae feedback affecting the slope 
of the predicted LF at the faint end. 

The above considerations suggest that the comparison 
of the observed and predicted UV luminosity density at 
z > 3 is critical at the bright end. Indeed, the details 



of the dust extinction are important at the bright end of 
the steep LF function since the dust extinction affects not 
only the observed UV luminosities but also the number 
of galaxies that are detected within the magnitude limit. 
Some spectroscopic information for the brightest sources 
can better constrain the dust content at these high z. 

For such a reason, the comparison between the predicted 
and observed UV luminosity densities should be performed 
over wider areas of the sky, to reduce fluctuations in the 
number of bright sources, and using different limiting mag- 
nitudes, extended to the faintest limits, to explore the dust 
content as a function of the galaxy star formation rate and 
mass. This will allow us to sample the shape of the LF 
at the faint end, which is crucial to assess whether com- 
plementary physical processes have to be included in the 
current hierarchical models. 



REFERENCES 



Arnouts, S., D'Odorico, S., Cristiani, S., Zaggia, S., Fontana, A., 

Giallongo, E., 1999, A&A, 341, 641 
Bertin, E., & Arnouts, S., 1996, A&AS, 117, 393 
Cole, S., Lacey, CO., Baugh, CM.; Prenk, C.S., 2000, MNRAS, 319, 

168 

Connolly, A. J., Szalay, A. S., Dickinson, M., SubbaRao, M. U., 

Brunner, R. J., 1997, ApJ, 486, Lll 
Efstathiou, G., Ellis, R. S., & Peterson, B. A., 1988, MNRAS, 232, 

431 

Ellis, R.S., CoUcss, M.M., Broadhurst, T.J., Heyl, J.S., Glazebrook, 

K. 1996, MNRAS, 280, 235 
Fernandez-Soto, A., Lanzetta, K. M., Yahil, A., 1999, ApJ, 513, 34 
Folkes, S., ct al., 1999, MNRAS, 308, 459 

Fontana, A., D'Odorico, S., Poli, F., Giallongo, E., Arnouts, A., 

Cristiani, S., Moorwood, A., Saracco, P., 2000, AJ, 120, 2206 
Fontana, A., Menci, N., D'Odorico, S., Giallongo, E., Poli, F., 

Cristiani, S., Moorwood, A., & Saracco, P., 1999, 310, L27 
Giallongo, E., D'Odorico, S., Fontana, A., Cristiani, S., Egami, E., 

Hu, E., McMahon, R. G. 1998, AJ, 115, 2169 
Heyl, J., CoUess, M., Ellis, R.S., Broadhurst, T. 1997, MNRAS, 285, 

613 

Lacey, C, & Cole, S., 1993, MNRAS, 262, 627 



Lilly, S., Trcssc, L., Hammer, F., Crampton, D., Le Fevre, O., 1995, 
ApJ, 455, 108 

Lin, H., Ycc, H.K.C., Carlberg, R.C., Morris, S.L., Sawicki, M., 

Patton, D.R., Wirth, G., Shepherd, C.W. 1999, ApJ, 518, 533 
Marzke, L.O. & da Costa, L.N. 1997 AJ, 113, 185 
Perlmutter, S. et al. 1999, ApJ, 517, 565 

Pettini, M., Kellogg, M., Steidel, C. C, Dickinson, M., Adelberger, 

K. L., Giavalisco, M., 1998, ApJ, 508, 539 
Poli, F., Giallongo, E., Menci, N., D'Odorico, S., & Fontana, A. 1999, 

ApJ, 527, 662 

Sandage, A., Tammann, G. A., & Yahil, A., 1979, ApJ, 232, 352 
Sawicki, M.J. Lin, H., Ycc, H.K.C. 1997, AJ, 113, 1 
Schmidt, M. 1968, ApJ, 151, 343 

Somerville, R. S., Primack, J. R., Faber, S. M. 2001, MNRAS, 320, 

504 

Steidel, C.C., Adelberger, K.L., Giavalisco, M., Dickinson, M. 

Pettini, M. 1999, ApJ, 519, 1 
Steidel, C. C, Giavalisco, M., Pettini, M., Dickinson, M., & 

Adelberger, K. L. 1996, ApJ, 462, L17 
White, D. A., & Fabian, A. C, 1995, MNRAS, 273, 72 
Zucca, E., et al. 1997, A&A,326, 477 



